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Abstract. A type of coupled variable coefficient modified Korteweg-de Vries system is 
derived from a two-layered fluid system. It is known that the formation, maintenance, and 
collapse of an atmospheric blocking are always related with large-scale weather or shorts 
term climate anomalies. One special analytical solution of the obtained system successfully 
features the evolution cycle of an atmospheric monopole type blocking event. In particular, 
our theoretical results captures a real monopole type blocking case happened during 19 Feb 
2008 to 26 Feb 2008 can be well described by our analytical solution. 



1. Introduction 

Atmospheric blocking is a fundamental large-scale weather phenomena in mid-high lat- 
itudes in the atmosphere that has a profound effect on local and regional climates. The 
life cycle of an atmospheric blocking always brings about large-scale weather or short terms 
climate anomalies. Therefore, it is rather important to predict an atmospheric blocking in 
regional midterm weather forecast and short-term climate trend prediction. 

There are three types of patterns for an atmospheric blocking anticyclone, i. e., monopole 
type blocking (or omega type blocking), dipole type blocking, and multi-pole type blocking. 
During the past many years, the dipole type blocking has been studied a lot since its first 
discovery by Rex [l| ■ Malguzzi and Malanotte-Rizzoli first used the Korteweg de- Vries (KdV) 
Rossby soliton theory to study dipole type blocking. While unfortunately their analytical 
results failed to describe the onset, developing, and decay of a blocking system. In fact, the 
important atmospheric blocking can be explained by many different theories except for KdV 
type equations. Recently, Luo et al. proposed the envelope Rossby soliton theory based 
on the deduced nonlinear Schrodinger (NLS) type equations and successfully explained the 
blocking life cycle numerically. More recently, we have found that variable coefficient KdV 
equation can analytically features the life cycle of a dipole blocking if introducing a time- 
dependent background field Q]. In addition, it has been revealed in Ref. that none zero 
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boundary values and time-dependent background westerly are vital factors to explain a life 
cycle of a dipole type blocking. 

In this paper, we are motivated to introduce time into the background flow and boundary 
conditions to derive new types of equations from a two-layered fluid to investigate atmospheric 
blocking systems. The paper is organized as follows. In Section 2, a type of coupled variable 
coefficient modified KdV type system is derived from a two-layered fluid model. Then we give 
a special analytical solution with many arbitrary functions and constants in Section 3. In 
Section 4, we assume some special values of the parameters in the analytical solution to obtain 
an approximate analytical expression for the stream functions which can describe a typical 
monopole type blocking event. The whole life cycle of the monopole blocking is graphically 
displayed, which really captures the feature of a real observational monopole blocking case 
happened in the 2008 snow storm in China. Last section is a short summary and discussion. 

2. Derivation of the coupled variable coefficient modified KdV type system 
The starting two-layered fluid model is 

qit + J{ih,qi} + pil>ix = o, (i) 

q2t + J{4>2,q 2 } + W2x = 0, (2) 

where 

Ql = Iplxx + tplyy + F(lp 2 - Ipl), (3) 

q2 = 1p2xx + Tp2yy + - 1p 2 ), (4) 

and J{a, b} = a x b y — a y b x . In Eqs. {I])-©, F is a weak coupling constant between two layers 
of the fluid and (3 = Po(L 2 /U), flo = 2(u>o/ao) cos(</>o), where ao is the earth's radius, ljq is 
the angular frequency of the earth's rotation and <j)Q is the latitude, U is the characteristic 
velocity scale. The derivation of the dimensionless equations (pQ) and ([2]) is based on the 
characteristic horizontal length scale L = 10 6 m and the characteristic horizontal velocity 
scale U = lOm/s Q|. A type of coupled KdV equations has been derived from the system 
CO)-®, and its Painleve property and soliton solutions have also been discussed 

aa. 

It is a common treatment to rewrite the stream function t/j by two parts, namely, ip = 
^oiVit) + V ; ' with i/jq being the background flow term, introduce the stretched variables 
£ = e(x — cot),r = e 3 t, (co is a constant), and take the background field ipo only as a linear 
function of y, and then expand the stream function as ij/ = YlnLi eTl ^i(£> Vi r )- Recently, 
considering the fact that background flow and the shear of the flow both could be time 
dependent, a new treatment was adopted in Refs. 0] to view background field rpo as a function 
arbitrarily depending on (y,t) and then further expand it as ipo = Uo(y) + Ylri=i en U n (y,T). 
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Below, a new type of coupled variable coefficient modified KdV type system is derived in 
a different way from that used in Ref. Q]. First, we introduce a variable transformation 
' — a\\x + a\2y,y' = a222/> and denote a| 2 = c\,a\ia22 = C2,a\ 1 + af 2 = C3. Then, 



x 



we rewrite the stream functions tpi = ipio(y', T ) + ^i(£> y' 1 T ){i = 1)2) with the stretched 
variables £ = e(x' — Cot),r = e 3 t, (co is a constant), and finally make the expansions ipxo = 
W,T)+En=i e "W,T),</>20 = f/o(y , ,r)+Er=ie n ^n(y , ,r), and $ = E~ 
In addition, we introduce F = Foe, (5 = A)^ 2 . Then we substitute all the expansions into Eqs. 
([1]) and ([2]) with Eqs. ([3]) and and then vanish all the coefficients of each order of e. For 
notation simplicity, the primes are dropped out in the following. 
In the first order of e, we obtain 

4>n = (C (T)t + A 1 (Z,T))B 1 (y,T) = (C Z + Ai)B 1 , (5) 
*hi = (C 1 (T)^ + A 2 ^,T))B 2 (y,T) = {C^ + A 2 )B 2 , (6) 



where B\ and B 2 satisfy 



(CO + C2V0y)Blyy ~ OlVtyyyyBx = 0, (7) 



and 



(CO + C2Uoy)B 2 yy ~ C 2 6 T QyyyB2 = 0, (8) 

respectively. The general solutions of Eqs. ([7]) and ([5]) read 

B x = (Fx J (c + cxV 0y )- 2 dy + F 2 ^ (c + cxV 0y ), 



and 



B 2 = \F 3 J (c + cxU 0y ) 2 dy + (c + ciC/ 0y ), 



(9) 



(10) 



where Fj, (i = 1, 2, 3, 4) are arbitrary integration functions of r. 
In the second order of e, we have 

^12 = Bs^ie + 5 4 A 2 + B 5 Al + fl 6 £Ai + B 7 Ai + B 8 £ 2 + B 9 £, (11) 

V>22 = #11^ + B12A1 + BxzA\ + B 14 ^2 + B 15 A 2 + £ 16 £ 2 + fli £, (12) 



where Bi = Bi(y,r), (i = 3,4, 16) are determined by a system of equations presented in 
Appendix A. 
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In the third order of e, if we assume 

^xa = B 19 A 2 ^ + B 2G A m + (5 2 iAi + B 22 £ + B 23 )A^ + B 2A A\ + B 25 Aj + S 30 ^^ 

+ (B 26 A 2 + £ 27 £ 2 + #28 + #29^1 + (B 31 £ + £32)^2 + + #34? + B 35 f 

+ J B 36 A% + (£37^2 + BnZAJAit + B 39 Al + (£ 40 £ + B 41 )A ± + B 42 A 2 d£, (13) 

^23 = -643^1^ + -E>44^2££ + (-B45^2 + -E>46<!; + -647)^2^ + -E>48^2 + B 4qA\ + B^A\ 
+ (£ 5 0^1 + B 51 f + B 52 + B 53 ^A 2 + (#55? + B m )A± + B 57 f + B 58 ^ + B 59 £ 2 

+ J B m A\ + (BetAx + B 62 ^A 2 )A 2 ^ + B 63 A 2 2 + (£ 64 £ + B 65 )A 2 + B m A 1( ^, (14) 

and requiring Bi = B; L (y,T),{i = 19, 20, 66) satisfy a system of equations (we do not 
write them down here for they are too long while can be easily retrieved following the above 
procedures), then we arrive at a coupled variable coefficient modified KdV type system 

Ai T + eie^i^ + (e i5 A\ + (e 42 ^ + e 24 )^i + e 32 £ 2 + e 30 + e 3 i£ + ei 2 A 2 )A^ + e 44 ,4 2 f 
+e 36 ,4 2 + (eu + e 15 £).Ai + e 21 A 2 ^ + (ds^l + e 2 o^2 + e 2 2^ + e 2 3)^2£ 

+(e 29 + en A 4 + e 28 £)Ai^ + e 38 £ + e 39 £ 2 + e 37 + ei 3 A 2 = 0, (15) 
A 2r + e 25 A 2 ^ + (e 35 A 1 + e 34 A 2 + (e 6 + e 7 £)A 2 + d + e 2 £ + e 3 £ 2 )^2£ + e 33 ^ 
+e 8 yl2 + (e 5 + e 4 £)vl2 + e 40 vli^ + (di£ + e 26 + d9^2 + e il A 1 )A 1 ^ 
+(e 9 £ + eio + e 27 ,4 2 )A 2? £ + e 43 ^i + e 47 £ + e 48 £ 2 + e 46 = 0, (16) 

with &i = ej(r), (i = 1, 2, ...,48) being arbitrary functions of the indicated variable. 

3. Special exact solutions 

Since equations (|15p -(|16 p constitute a coupled variable coefficient nonlinear system, it is not 
easy to obtain its general solution. Here we present a quite special solution of the system. 
It is easy to see that if we suppose A\ = aA 2 with constant a, then Eqs. (|15p and (|16p 
degenerate to one modified KdV type equation 

A 2r + e w A 2 ^ + (mio + ae 17 A 2 + e 2sO A 2^ + ae^A^ + ae 3 §A\ + (™5 + e 4 £)^2 
+(a 2 e 45 ^2 + (ae 42 £ + m 6 )A 2 + e 32 <!; 2 + m 2 £ + mi)A 2? + e 48 £ 2 + e 47 £ + e 46 = 0, (17) 

where rrn = mj(r), (i = 1, 2, 5, 6, 10) are given by mi = e\ + ae 2 Q,m 2 = e 2 + aen, ms = es + 
ae 43 , m 6 = e 6 + ae 35 + a 2 e 4i + aei 9 , mi = do + ae 40, when e 3 = e 32 , e 9 = e 28 , ei 5 = e 4 , e 34 = 
a 2 e 4 5,e 33 = ae 44 ,e 7 = ae 42 ,e 37 = ae 46 ,e 38 = ae 47 ,e 2 7 = ad7,e 39 = ae 48 ,e 8 = ae 36 ,e 25 = 
d6,e 2 i = ado + fl2e 40 - ae 29 ,d 3 = -aeu + ae 5 + a 2 e 43 ,e 2 3 = -ae 30 + ae% + a 2 e 2 6,e 2 2 = 
a 2 en — ae 3 i + ae2, e2o = — «d8 — aei2 — a 2 e2 4 + ae6 + a 2 d 9 + a 2 e35 + a 3 e 4 i- 
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Further, it can be verified that the modified KdV equation (fT7|) can be transformed to the 
standard one 



P T + 6P 2 P X + Pxxx = 0. (18) 

if 

A 2 = -^p/l + f 2 P(X,T), (19) 
where T = / 4 (r), X is given by 

and fi = fi(r), (i = 1,2,3,4) are arbitrary functions, with conditions 

mlfi + e i5 f 3T fi T fa a 6 f iT 

m 1 = " ,m 2 = — r ,m 5 = — — , ei 6 = — 777 (21) 



e 45 = — g , 2a e 45 e 46 - m 5 m 6 e 45 - m 6r e 45 + m 6 e 45T = 0, (22) 
/1/2 

ei7 = e 2 8 = e 32 = e 36 = e 4 = e 42 = e 44 = e 47 = e 48 = mi = 0. (23) 



Hence, it is easy to obtain exact solutions of Eq. (|T7j) based on the solutions of Eq. (|T8[) 
through Eq. (fl~9|) with Eq. ([20]) , One classical soliton solution of the mKdV equation (fl~8|) 
is P = _ffsech(l^X — K 3 T) with a constant -RT. Now just using this typical solution, we can 
easily write down a special solution of the original system {I])-© as 



ipi « V (a 2 2y,T) + eVi(a 22 y,r) + eB 1 (a 2 2y,r) |C e(anx + a i2 y - c t) 
+a/ 2 Asech ^ 2 (anx + a i2 y - c t) — K°f 4 

%Ij 2 w U (a 22 y, t) + eUi(a 22 y, r) + eB 2 (a 22 y, r) |de(aiix + a i2 y - c i) 

, [6-fC/ie 3if/ 3 3 ' 
+/ 2 Asech ^ 2q2 (anx + a i2 y - c i) ^4 K 



m G fif 



2af 



At 



(24) 

2a?f 



At 



(25) 



with r = e 3 t, I?i and B 2 determined by Eqs. ([9]) and (I10p . respectively. From the analytical 
solution of stream functions ipi and -02; it is easy to derive the corresponding background 
westerly flows u\ = —di^i^/dy = —dVo(y,T)/dy — edVi(y,T)/dy, and u 2 = —dip 2 o/dy = 
—dUo(y,T)/dy — edUi(y,r)/dy, which are all left as arbitrary functions of the indicated 
variables. 
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4. Analytical diagnosis 



By selecting the arbitrary functions and constants appropriately, the approximate analyti- 
cal solution (|24 p - (|25p can be responsible for different kinds of atmospheric blocking phenom- 
enon. Here we give a typical example when the functions and constants are taken as 

V = -0.006y 2 + O.OOly - 3sech(0.25y - 1.25), a = C = h = f 2 = a 22 = 1, 

c o = m 6 = F 1 = V 1 = 0, f 3 = -24/ 4 (r) + 0.1, F 2 = sech(0.5(i - 4.6)), (26) 

a n = 0.3, ai2 = -0.06, c x = 10, e = 0.1, K = -6. 



In this case, the solution (|24j) can describe a monopole type blocking event as depicted in 
Figure [TJ Considering the fact that the basic flow, mainly the basic westerly and the shear 
of the basic westerly, plays a significant role on the blocking developing process l|, 3, S, Q], 
it is hence reasonable to introduce the y 2 related term. 



(a) Day 1 



(b) Day 2 



(c) Day 3 




Figure 1. A monopole type blocking life cycle from the theoretical solution 
with the parameters given by (|26p. The contour interval (CI)=0.12. 



Evidently, a whole life cycle of a monopole type blocking event, namely, the onset, develop- 
ment, maintenance, and decay processes, are clearly presented in Figure [TJ The streamlines 
are gradually deformed, and the anticyclonic high in the north develops at the second day 
(Fig. la). It is strengthened daily. At around the fourth day (Fig. le), it is at its strongest 
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stage and then become weaker and eventually disappear after the seventh day (Figs. le-h). 
Obviously, Fig. 1 possesses the phenomenon's salient features including their spatial-scale and 
structure, amplitude, life cycle, and duration. Therefore, Fig. 1 is a very typical monopole 
type blocking episode in one layer of the fluid. 

A real observational blocking case happened during 19 Feb 2008 to 26 Feb 2008 is shown 
in Fig. [5] appearing a monopole pattern. It is easily found that the life cycle of this blocking 
lasts about eight days experiencing three stages: onset (19-20 Feb, 2008), mature (21-22 Feb, 
2008) and decay (23-26 Feb, 2008) periods, which resembles those displayed in Figure 1. 



(a) 19 Feb. 2008 (b) 20 Feb. 2008 (c) 21 Feb. 2008 




(d) 22 Feb. 2008 (e) 23 Feb. 2008 (f) 24 Feb. 2008 




(g) 25 Feb. 2008 (h) 26 Feb. 2008 




Figure 2. Filtered geopotential height at 300-hPa pressure level of a mono- 
pole blocking case during 19-26 Feb 2008. The x-axis is longitude, and the 
y-axis is latitude. Contour interval is 6 gpdm. 

5. Summary and discussion 

Considering a time-dependent basic westerly, we have derived a type of coupled variable 
coefficient modified KdV type system from a two-layered fluid model, with the time depen- 
dent coefficients resulted from the time dependent basic flows and time dependent boundary 
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conditions. Instead of possessing a linear meridional shear, the mean flow is assumed to be 
a combination of a quadratic function of y and sech(y). The boundary conditions remain as 
unknown functions with some complicated relations. Under a set of parameters, our analyt- 
ical solution nicely described a real observational blocking life cycle from 19 Feb to 26 Feb 
2008, indicating the onset, mature and decay phases in its developing process in one of the 
layers of the fluid. In the other layer, similar or different types of blocking can be found when 
setting the unknown parameters. 

Therefore, it is revealed that blocking can also be governed by modified KdV type system 
and its analytical solutions can also features the life cycle of a blocking. It is worth further 
investigations on how variations of the time-dependent background westerlies influence the 
type of a blocking and the evolution of a blocking during its life period since the role of the 
weak westerlies on a blocking is still unclear. 
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Appendix A 

Bi = Bi(y, t), (i = 3, 4, 18) are determined by the following equations 

c\(cq + C2V 0y )B 3y y - C2a\ 2 Vo y yy B ' 3 + 2a yi (ci a 1 1 V 0y + c a 2 2)Biy = 0, (27) 
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ci(c + c 2 Voy)B iy y - cic 2 VoyyyB 4 + (c + c 2 V 0y )F B 2 = 0, (28) 

2(ciclV^ y + 3c cic 2 V 2 y + CqCi + 3(%c 1 c 2 V 0y )B 5yy + ci(^BlV Gyyy V Gyy 

-2cic 2 (2coc 2 V y + eg + (^Vo y )V Qyyy B 5 - c 2 (c + cic 2 V Qy )VQ yyyy Bf = 0, (29) 

( c i c 2^o 3 j/ + Scocic^oy + CqCi + Sc^clclVo^Beyy - cic|(c C + c^oV^Vc^^ 2 
+cic 2 c 2 BfCoVo yyy Vo yy — cic 2 {clV^ y + ZeoCzVoy + cl)V 0yy yB 6 = 0, (30) 

(^o 2 j/ a ii a 22 + 2ciVbj,a2 2 anco + c\c^)B 7yy - cjc 2 (c + c 2 V y)Vo y yyB 7 + a^a^BiVoyyyViy 
-ci(cq + coc 2 U 0y + c c 2 V y + clt/oyV^-Fo-Bi - cfc 2 (co + c 2 Voj / )Vi W j,.Bi = 0, (31) 

2c?(ciy 3 y + 3co^Vg, + eg + 3^02^)^ - 2c 2 c 2 (c 2 1 + 2c c 2 y 0y + c^^o^s 

— CiC 2 Cq(CQ + C 2 VQy)V0yyyyBf + c\c\B\ C^V§yyyV$yy = 0, (32) 

(VQycjcl + 2ciy j/Cocic 2 + Ci(^)5io w - cjc 2 Co(co + QVb^Vi^Bi - 2coF B 2 C 1 cic 2 V 0y 
-c 2 c 2 (c -B 9 + c 2 B 9 V 0y - c 2 C BiVi y )Voyyy - C cic 2 (co + c 2 Vb y )iWoy-Bi 
+ciCi(cg + c^o^FoSa - c ciCo(c + ca^o^o^i - c?(c + caV&^Vfowr = 0, (33) 

ci(c + c 2 Uo y )B 1 iyy + 2a i2 (d ai i [/ j, + a 22 co)B 2y - cic 2 U 0yyy B u = 0, (34) 

ci(c + c 2 U 0y )B 12 - anai 2 ^0j/ ra -Bi2 + (co + c 2 C/ j/)^b^i = 0, (35) 

2c 2 + 3co4U$ y + 3cgc 2 ?7os, + 4)B 13yy + +44 B 2 U Oyyy u Oyy 

-2c\c 2 {clUQy + CQ + 2C0C 2 U0y)B 13 U Q yyy + C\C\(CQ - C 2 U y)B 2 U0yyyy = 0, (36) 

CiO^oj/ + 3c oc|^oj/ + 3cgc 2 f7oj, + c^-Buyj/ - c?c|Ci(co + c 2 U 0y ) B^U 0y yyy 

-<?\C 2 ((? 2 XJly + 2c c 2 [/ j/ + cl)B 14 Uoyyy + c 2 c 2 Ci ?7os/j/s/ ^oyy = 0, (37) 

c i( c 2^% + 2c 2 £/os/C + (^)Bi5 J ,j / - cic 2 (c + c 2 U 0y )B 2 Uiyyy - c 2 (c + c 2 [/o2/)-FoVoj/-B 2 

-C1C 2 (C + C 2 U y)U0yyyBl 5 ~ Cq(CQ + C 2 U y)F B 2 - ClclB 2 U0yyyUly = 0, (38) 
2C l( U 0y C 2 + SU^yClcO + 3U yC 2 cl + cl)B 17 yy - clclCf(C0C 2 U y)B 2 2 U0yyyy 

+cjcjB 2 CfUoyyyUoyy - 2c\c 2 (2c 2 U()yCQ + CQ + clUQy)U yyyB w = 0, (39) 

c i( c 2 U Q y + 2coc 2 [/oj/ + cl)Bi 8 yy - ci(c + c 2 U y)Uoyy T - c\c 2 Ci(c + C 2 U y)B 2 Ulyyy 
+cic 2 (c + c 2 U y)UoyyyB lo + c^Fq{CqB 1 - CiB 2 ) + c 2 Co(c 2 [/ j/ + 2c )F 5i?7oj / 
-c 2 Ci(c + c 2 V y)F U yB 2 - co^FoB^oyC 1 ! + cx(? 2 C x B 2 UQ yyy U Xy = 0. (40) 
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